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Atomic-resolution scanning tunneling spectroscopy is used to study the spectrum of quasiparticle 
states in superconducting SnRuCU. The measured temperature dependence of the quasiparticle 
spectrum is consistent with a multi-band superconducting order parameter with one band exhibiting 
a line of nodes, if a constant density of states offset can be ascribed to a surface effect. In magnetic 
fields H < H C 2, a vortex lattice, square and oriented parallel to the (110) direction of the Ru02 
lattice, is detected by spectroscopic imaging. Each vortex exhibits one flux quantum and a strong 
zero-bias conductance peak. 

PACS numbers: 68.37.Ef, 74.20.Rp, 74.70.Pq 



Sr2RuC>4 displays an exotic form of superconductivity 
at temperatures below 1.45 K 0,0. Because of the inde- 
pendence of the spin susceptibility on temperature below 
T c QL , the Cooper pairs are believed to be in a spin 
triplet 'equal spin pairing' (ESP) state. Muon spin rota- 
tion studies also show a time-reversal-symmetry breaking 
(spontaneous magnetization) signal which is also consis- 
tent with spin-triplet superconductivity 0]. 

Many discussions have proceeded in terms of a purely 
2-dimensional version of the p-wave order parameter 
(OP) d = Az(k x ± iky). But it now appears that the 
Sr2Ru04 OP is somewhat more complex than this. It 
has recently emerged that (i) the low temperature elec- 
tronic heat capacity C e , (ii) the nuclear spin relax- 
ation time 0, (iii) the thermal conductivity 0, (iv) the 
London penetration depth [lOj and (v) the ultrasound 
attenuation |ll| have low temperature power laws which 
are consistent with the presence of line nodes (Fig. |TJ)) . 

A microscopic model consistent with these observa- 
tions has also recently emerged. In Sr2Ru04, three elec- 
tronic bands cross the Fermi level leading to three closed 
Fermi-sheets in momentum space: a (X centered hole- 
like), (3 and 7 (r centered electron like) 0, 0|. It is 
proposed that on the 7 sheet the superconductivity is 
dominant while the remai ning two sheets have induced 
superconducting gaps 0, Il5| . Although this model of 
Sr 2 Ru04 is elegant, plausible, and consistent with ex- 
periments, it is by no means certain |2(. 

Tunneling spectroscopy has historically played a key 
role in the study of superconductivity [lj| . This is be- 
cause the differential tunneling conductance g — dl/dV 
at bias voltage Vb between a sample and metal tip is 
proportional to the sample local density of quasipar- 
ticle states (LDOS) at energy E = eVb] g(r,Vb) oc 
LDOS(r,i? = eVb). Scanning tunneling spectroscopic 
imaging has further advanced the utility of this technique 
by allowing mapping throughout space of <?(r, Vb). For 
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FIG. 1: a) Crystal structure of Sr2Ru04. b) Schematic of 
various gap order: top is a p-wave, middle is vertical line 
nodes and at the bottom is a multi-band model with both a 
modulated p-wave and vertical line nodes on different Fermi 
surfaces. 



example, in s-wave NbSe2, the vortex lattice and core 
structure have been determined with this technique |l7| . 
and in d-wave cuprate superconductors it has revealed 
the vortex structure 0], nanoscale electronic disorder 
[Tij. impurity atom effects [2(|, and quasiparticle quan- 
tum interference |2l| . In this paper we apply these tech- 
niques to Sr2Ru04. 

Several important STM studies of Si'2Ru04 have al- 
ready been reported. At temperatures high compared 
to T c , topographic imaging of surface structure reveals 
unexpected reconstructions of SrO surface j^- Tunnel- 
ing spectroscopy has been carried out below T c but on 
an unidentified, and ramified, surface yielding the unex- 
pected result that the LDOS at E = was about 90% 
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FIG. 2: a) Topographic image of a Ti doped sample showing 
the usual square and uniform lattice (and some Ti impuri- 
ties), b) Spectra taken on the surface a) away from a Ti 
atom showing a 50 mV and a 10 mV gaps. The inset shows 
that there is no superconducting feature seen at low energy, 
c) Topographic image taken on the pure sample studied here 
showing a very different surface structure, d) spectra taken 
on the image of c) and showing a very different V-like spectra 
with weak features (small slope changes) at 50 and 10 mV 
and reaches a very small conductance around zero bias. The 
inset shows data taken at small bias voltage (±1 mV) showing 
the superconducting gap. The topographs were taken with a 
setpoint of 100 mV and 0.05 nA. 



of normal value even at low temperature |23j . A possible 
explanation was the level of disorder seen in the images. 
Nonetheless, this study was an important advance since 
it detected a partially gapped quasiparticle spectrum. 

We use a dilution-refrigerator based STM which 
reaches 15 mK in fields up to 8.5 T with atomic- 
resolution scanning tunneling spectroscopy. This sys- 
tem also allows samples to be cleaved in cryogenic ul- 
tra high vacuum inside the vacuum space of the dilu- 
tion refrigerator and then immediately inserted into the 
STM head mounted below the mixing chamber. We use 
floating-zone-fabricated Sr2Ru04 single crystals which 
have T c =1.45 K. 

In Sr2Ru04, this cleavage technique can, in principle, 
result in two c-axis normal crystalline surface layers: the 
SrO layer and the RuC>2 layer (see crystal structure in 
Fig. I^i. Topographic images of two surface types we ob- 
serve are shown in Fig. and Fig. The first surface 
layer was identified as being SrO using Ti atom substitu- 
tion experiments |2 1| . This surface never exhibits a su- 
perconducting quasiparticle spectrum but instead shows 
the spectrum in Fig. 12b 

A topographic image of a second type of cryo-cleave 



surface (likely RUO2) is shown in Fig. This image 
shows a quasi-one-dimensional structure on this surface. 
The width of the "bands" in the figure are 4ao (where ao 
is the in-plane lattice constant, 0.386 nm) while the peri- 
odicity of the "beads" running the length of the "bands" 
is 2ao- Domains of both orientations of the quasi 1-d to- 
pographic structure are clearly seen everywhere on this 
surface. Surface reconstruction effects that are sensitive 
to the cleavin g re cipe have been seen in some measure- 
ments |U I25I l2(i| . It is not clear how these high tem- 
perature cleaved reconstructions are related to the one 
reported here. 

Independent of structural details of this RuO layer 
reconstruction, the surface shown in Fig. [2t displays a 
differential tunneling conductance spectrum as shown in 
Fig. [2Ji, and with high resolution in Fig. [3^). We have 
studied the spectrum over regions up to 1.1 /1m 2 and it 
is highly repeatable and reliable, having almost identi- 
cal characteristics independent of exact atomic details of 
surface reconstruction. 

We carried out the following tests to demonstrate that 
these spectra represent the superconducting quasiparti- 
cle spectrum of S^RuO^ First, the temperature de- 
pendence studies show the observed gap disappears at 
T c (see Fig. Et>) filling from 40% to 100% of normal 
LDOS(i? = 0). Second, at millikelvin temperatures, this 
gapped spectrum has disappeared at a magnetic field of 
/io-ff = 100 mT, which is very close to the c-axis H C 2 
of the superconducting state. Third, in the presence of 
fields H < H C 2 this spectrum shows the characteristic 
spatial variations expected from the core of vortices ex- 
hibiting 0o of magnetic flux (see Fig. . We therefore 
conclude that we are dealing primarily with the quasi- 
particle spectrum of superconducting OP of Sr2Ru04. 

To analyze the shape of the spectra and their temper- 
ature dependence we consider three different models of 
OP; an isotropic p-wave gap 



A(0) = Ag, 



a gap with a vertical line of nodes 
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and a modulated p-wave gap (which is combined with 
vertical line nodes for the multi-band model): 



A(</>) 



a p-mod ±%4 
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In this latter case, the in-plane modulation of the gap 
along cj> has peak-to-peak amplitude of a (i.e. a maximum 
of 1 and a minimum of 1 — a). 

The local conductance g(r, V&, T) is obtained from the 
local density of states of the various gap models (Eq. la- 
ic) by including thermal smearing T. Therefore predic- 
tions for the differential tunnel conductance of the gaps 
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FIG. 3: a) Shows three different fits to the lowest temperature 
conductance data. The dotted (blue) line is for an isotropic p- 
wave gap. The dashed (red) line is for vertical line nodes and 
the solid (black) line is for the multi-band model discussed 
in the text, b) Temperature dependence of the data and the 
associated multi-band model fits. The different temperatures 
are shifted vertically for clarity, c) The A p _ mo( j gap param- 
eter obtained from the fit in b) compared with a BCS like 
temperature dependence shown by the line (Aii nc is kept a 
constant fraction of A p _ mo d hence has the same temperature 
dependence) . 



are obtained as g P (V) (which could as well represent an s- 
wave OP), giinc(V) and <7 p - mo d(^0- Note that horizontal 
line nodes (An cos(k z c)) gives the same conductance as 
vertical line nodes, and that the orientation of the nodes 
as well as the exp(±j(/)) part of the gaps also does not 
modify the conductance. 

To relate the three models to the measured conduc- 
tance g(r, Vb) in the experiment (Fig.|3i), we fit the data 
to a general form which includes a background density of 
states which is added as a fraction / of the total conduc- 
tance, 



g(V)=A (l-f)G(V) + f 



(2) 



and for the multi-band model 



9(V) =A{(1- /) [(1 - X) 5 Hne(T0 + a^mod(V0] + /} 

(3) 

where V is the bias voltage, A is the amplitude of the 
conductance and G is either gn nc or g p . We consider /, 
which is about 40%, to be a spurious effect. This is rea- 
sonable since heat capacity experiments show the number 
of thermally occupied states in bulk falling precipitously 
at these low temperatures. In case of the multi-band 
model, x is the fraction of the non-background LDOS 
due to the modulated p-wave gap (the 7 band). There- 
fore we have 2 or 3 adjustable parameters (A, f and x) 
in addition to the predicted g(V) of the gap function of 
interest. The temperature is fixed to the measured sam- 
ple temperature, except at very low temperatures below 
250 mK where we fix the temperature to 250 mK be- 
cause (1) we believe from independent studies that this 
tip did not cool below that temperature and (2) this is 
the sample temperature at which the measured spectra 
start changing. 

As can be seen in Fig. the simplest p-wave fit 
G = g p is quite poor compared to the line nodes. The 
multi-band model is better than a pure line nodes model 
(possibly because of the larger number of parameters). 
Other gap-models with extra parameters over a simple 
line node also fit as well, such as a modulated horizontal 
line of nodes: A cos(fc z c) [1 — /?cos(4<^)] (but this would 
fail to reproduce the heat capacity jump at T c ). The 
particular multi-band model was selected because it is 
consistent with a recent study of the dependence of the 
heat capacity on the angle of the magnetic field [27j . We 
caution that the value of x is controlled not only by the 
density of states of the different bands but also by the 
tunneling matrix element. 

We next study the temperature dependence of these 
tunneling spectra. In Fig. [SJd, we show averaged spectra 
taken at the 17 temperatures between 20 mK and 1.6 K. 
They were obtained after averaging about 100 individual 
spectra (25 s each) with a bias modulation of 20 //Vrms 
and setpoints of 10 mV / 0.25 nA or 5 mV / 0.15 nA. 
Different locations were used for data acquisition at each 
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temperature. The solid lines in Fig. |3Jd are the result of 
fits to the multi-band model. The lowest temperature is 
fitted the same way as in Fig. [3]} above. For all other 
temperatures /, x and a are kept constant at their low 
temperature value (0.38, 0.30 and 0.37 respectively), T 
is set to the measurement temperature but Ag~ mod is left 
to evolve while keeping the ratio Ag~ mod /A ino fixed to 
the low temperature value of 1.9. 

From all these fits, a temperature dependence of the 
gap AQ _mod can be extracted and the result is shown in 
Fig. |2t where the solid curve is a fit to a BCS temperature 
dependence of the gap. A very similar graph is obtained 
if a vertical line nodes model is used. 

A second important use of spectroscopic mapping is in 
the determination of effects of magnetic field on the su- 
perconducting OP. This can be very revealing since the 
structure of the vortex lattice and the electronic struc- 
ture of the core are closely related to the symmetry of the 
OP. For the case of as-wave OP, a number of theoret- 
ical studies [2^, |2!J, |3fJ, |31| have explored the predicted 
structure of the vortex core in 2-d OP d = Az(k x ±ik y ). 

Experimentally, the flux lattice has been measured by 
small angle neutron scattering J^]. However, until the 
present work, nothing is known about the local electronic 
structure of Sr2Ru04 vortex cores. To explore these is- 
sues we mapped point spectra at >140 individual loca- 
tions in B < fioH C 2- We detected 5 vortices with mean 
separation between cores of 180±30 nm. Within our sta- 
tistical uncertainty, the lattice is square and oriented par- 
allel to the (110) axes of the crystal lattice which is quite 
consistent with previous neutron studies. We detect the 
vortices from the changes in g(r, 0) at the vortex center, 
forming a pattern consistent with a square lattice and 
the separation of the core locations at B = 50 mT yields 
a flux $ = 1 ± 0.15<I>o threading each core. 

To study the Sr2Ru04 vortex core itself, we measure 
g(r, V) at regions nearby the core (Fig. ^Ji) and along 
lines connecting two cores (crystallographic (110) direc- 
tion) as shown in Fig. 2t>. The most striking feature is 
the very strong zero bias conductance peak (ZBCP) of 
FWHM 250 ^eV and maximum intensity 115% of the 
normal density of states. The coherence length can also 
be estimated from decay length of the gap suppression 
away form the core to be 50 nm in reasonable agreement 
with the estimate of 65 nm from penetration depth stud- 
ies. 

A number of important conclusions can be drawn. 
First, the RuO surface of Sr2Ru04 exhibits a gapped, 
particle-hole symmetric, density of states spectrum con- 
sistent with a superconductor. Second, the temperature 
dependence of the observed g(r, E) is best fit by a super- 
conducting multi-band gap with at least one band sup- 
porting a line of nodes [3Jj . This provides further corrob- 
oration for this picture of the gap structure in Sr2Ru04. 
Third, direct imaging of the vortex lattice by spectro- 
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FIG. 4: a) Averaged point spectra at the center (red) and in 
between (black) vortex showing the ZBCP at the center, b) 
Line cut through a vortex (as shown in the inset). The inset 
also shows the location of the 5 places where we identified a 
vortex core. The data was taken with a setpoint of 2 mV and 
0.15 nA. 



scopic mapping is now possible. Because of this, the long 
mooted experiment to search for chiral domains that can 
exist in a k x ± ik y OP [U becomes technically feasible. 
Finally, the large ZBCP discovered in Sr2Ru04 is consis- 
tent with p-wave models that did not consider line nodes. 
It will require new theoretical models to understand the 
ZBCP in a multi band model. 
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